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Abstract: In some areas of cultivation, a lack of salt tolerance severely affects plant 
productivity. Apple, Malus x domestica Borkh., is sensitive to salt, and, as a perennial 
woody plant the mechanism of salt stress adaption will be different from that of annual 
herbal model plants, such as Arabidopsis. Malus zumi is a salt tolerant apple rootstock, 
which survives high salinity (up to 0.6% NaCl). To examine the mechanism underlying this 
tolerance, a genome-wide expression analysis was performed, using a cDNA library 
constructed from salt-treated seedlings of Malus zumi. A total of 15,000 cDNA clones were 
selected for microarray analysis. In total a group of 576 cDNAs, of which expression 
changed more than four-fold, were sequenced and 18 genes were selected to verify their 
expression pattern under salt stress by semi-quantitative RT-PCR. Our genome-wide 
expression analysis resulted in the isolation of 50 novel Malus genes and the elucidation of a 
new apple-specific mechanism of salt tolerance, including the stabilization of photosynthesis 
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under stress, involvement of phenolic compounds, and sorbitol in ROS scavenging and 
osmoprotection. The promoter regions of 111 genes were analyzed by PlantCARE, 
suggesting an intensive cross-talking of abiotic stress in Malus zumi. An interaction network 
of salt responsive genes was constructed and molecular regulatory pathways of apple were 
deduced. Our research will contribute to gene function analysis and further the 
understanding of salt- tolerance mechanisms in fruit trees. 

Key words: salt stress; Malus zumi; microarray; salt-responsive gene; interaction network 



1. Introduction 

Salt stress severely affects crop production worldwide. Efforts have been made to uncover the 
mechanisms of salt tolerance in model plants using molecular and genomic approaches. Many salt tolerant 
genes have been identified [1,2], which release plants from salt stress through stress perception, signal 
transduction and transcriptional regulation for cellular responses including ROS scavenging, osmolyte 
accumulation, and transportation of water and ions through the plasma membrane and tonoplast [3,4]. The 
over-expression of a number of salt-induced genes confers stress tolerance to the transgenic plants [5,6]. 

Three primary signaling pathways in salt-stress response (CDPK, CIPK/SOS, and MAPK pathways) 

2_|_ 

have been reported in plants [4]. The Ca -dependent signaling (CDPK pathway) is extensively studied, 
which induces the expression of DRE/CRT transcription factors and other types of LEA-like genes. The 
CIPK/SOS pathway appears to be relatively specific for the ionic transportation, in which high Na + 
stress initiates a calcium signal that stimulates the SOS3-SOS2 protein kinase complex, which then 
activates the Na + /H + exchange transporter SOS1 and regulates the expression of several other 
salt-responsive genes. In addition, SOS3-SOS2 may activate or suppress the activities of other 
transporters involved in Na + homeostasis. The MAPK pathway regulates production of compatible 
osmolytes and antioxidants, and may also participate in cell cycle regulation under salt stress. 

Although the molecular mechanism of salt tolerance has been intensively studied in model plants, in 
perennial woody species such as apple, the mechanism remains unclear [7,8]. Apple {Malus x domestica 
Borkh.) is one of the most valuable horticultural fruit crops in the world. While it produces the highest 
yield in China, it is subject to severe salt stress in many areas [9,10]. A better understanding of the 
genetic bases of salt tolerance would contribute to the molecular breeding of salt resistant apples. 
Rootstocks play a key role in salt tolerance. The increasing secondary salinity threats apple production. 
Because in the commercial orchards, all the apple cultivars are grated onto rootstocks, which are 
subjected to the high salinity in the soil. Malus zumi is a local rootstock, distributed in the Northeast of 
China, which can survive high salinity (approximately 0.6% NaCl). 

The draft apple genome has been recently published, with the majority of the genes not functionally 
annotated [11]. Furthermore, apple transcriptomics have been performed using microarray analysis [12], 
however no analysis or isolation has been made of ESTs from samples under salt stress. To gain a 
genome -wide view of salt response, a cDNA library was constructed in Malus zumi under salt stress by 
SMART . A total of 15,000 cDNA clones were selected randomly for microarray analysis, among 
which, 576 cDNA clones changing their expression significantly (more than four-fold) under stress were 
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sequenced. Our research resulted in novel gene isolation and new molecular insights of apple-specific 
salt-tolerant mechanisms in term of ROS scavenging, osmoprotection, and photosynthesis in Malus zumi. 
The czs-elements in the promoter region of salt-responsive genes were analyzed. A network of salt 
responsive genes was constructed and molecular regulatory pathways of apple were deduced. This 
helped to elucidate salt-tolerant mechanisms from a novel perspective. Our research will contribute to 
identification of genes in response to salt-tolerance and as well as providing new insight into 
understanding salt-tolerance mechanisms in apple trees. 



2.1. Salt-Responsive Genes Identified by Microarray Experiment and Semi-Quantitative 
RT-PCR Validation 

According to three screening criteria, a total of 2562 cDNAs were isolated, whose expression 
changed more than two-fold under stress (Figure 1). Their corresponding genes were defined as 
salt-responsive. This suggested that approximately 17% of the tested cDNAs were salt-responsive, 
similar to results reported in wheat [13]. A total of 1713 cDNAs were up-regulated and 849 cDNAs were 
down-regulated, among which the expression of 1431 cDNAs changed two- to three-fold, 493 cDNAs 
changed three- to four-fold, 245 cDNAs changed four- to five-fold, and 393 cDNAs changed greater 
than five-fold (Figure 1). There are accumulated reports that genes which observe a fold-change of two 
may show no, or very little, differential expression due to the technical problems. The genes which 
showed changes in expression of over four-fold might be more possible to be involved in salt stress 
response in Malus zumi. As Zhao et al [14] describe, The 235 cDNAs which changed four- to five-fold, we 
functionally classified, according to the Gene Ontology (GO) prediction, into 1 1 catagories (Table 1): 
signal transduction (6%), ROS elimination (6%), osmoprotection (5%), cell maintenance and 
development (6%), photosynthesis (16%), transporter (5%), metabolism (17%), stress tolerance (7%), 
protein related (7%), others (4%), and unknown (21%) (Table 1). A complete list of the functional 
classification was shown in Supplementary file (Table SI). 

Figure 1. Distribution of genes according to the expression change in response to salt treatment. 
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Table 1. The functional categorization of the putative salt-responsive genes. 



Functional category 



Percentage of unigenes (%) 



Photosynthesis 
Transporter 
Metabolism 
Stress tolerance 
Signal transduction 
Protein related 
ROS elimination 
Osmoprotection 
Cell maintenance and development 
Others 
Unknown protein 



21 



17 



16 



4 



7 



7 



6 



6 



5 



6 



5 



The microarray results were validated by semi-quantitative RT-PCRs performed on 1 8 selected genes 
(Figure 2). MzTSL, MzILR, MzCIPK6, MzIAA, MzSCL, MzSTO, MzGTL, and MzDREBl were regulatory 
genes that function in the upstream. MzRCA, MzLhcb2, MzRPE, MzSDH, MzHSP70, MzNTR, MzRD22, 
MzPIP2, and MzTIP2 were functional genes at the downstream level. As the maintaining photosynthesis 
efficiency under salt stress might be the one of the major mechanisms in Malus zumi, we chose three 
genes {MzRCA, MzLhcb2, MzRPE) involved in photosynthesis. In addition, we also select an unknown 
gene MzUKl for expression confirmation. In the microarray, MzTSL, MzILR, MzIAA, MzSTO, MzRCA, 
MzLhcb2, MzHSP70, MzNTR, MzRD22, MzPIP2, MzTIP2, and MzUKl were induced by salt stress, 
when MzDREBl, MzGTL, MzSCL, MzCIPK6, and MzSDH were suppressed by salt stress 
(Supplementary file (Table SI)). With the exception of MzNTR (sodium dependent phosphate 
transporter; Bankltl498107), MzDREBl (AP2 transcription factor; Bankltl494547), MzCIPK6 (CIPK6; 
Bankltl494666), and MzGTL (GT-like trihelix DNA-binding protein; Bankltl498102), the RT-PCR 
results of the other 14 tested genes (78%) were consistent with the microarray data. Similar levels of 
validation are reported in other microarray studies [15]. This demonstrated the satisfactory quality of our 
experimental procedures. 

Among the eight genes encoding regulatory proteins, MzSCL (GRAS family transcription factor; 
Bankltl494667) was down-regulated and the other seven genes were up-regulated by salt-stress. MzSCL 
belongs to GRAS transcription factor family. Its expression is reduced steadily under salt stress (Figure 2). 
The seven salt-inducible genes could be divided into two groups according to the time of peak 
expression. The first group included MzSTO (SALT TOLERANCE Homolog Protein; Bankltl493633), 
MzDREBl, MzILR (IAA-LEUCTNE RESISTANT3; Bankltl495536), and MzTSL (Tousled-like 
Serine/threonine Kinase; Banklt 1494726), which responded to salt stress immediately after treatment 
and reached the peak expression as early as two hours (Figure 2). MzSTO, MzIAA (Bankltl493693), and 
MzGTL (Bankltl498102) of the second group were up-regulated after treatment and reached the 
expression peak at four hours (Figure 3). The early response of the tested transcription factors and 
kinases suggested their possible roles in salt-stress signaling. 
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Figure 2. The validation of the expression of 18 selected salt-responsive genes under salt 
stress by semi-quantitative RT-PCR. The expression levels of tested genes were normalized 
according to the corresponding actin amplifications, and were presented under the bands. 
The peak of expression level was colored with red. Except for MzNTR, MzDREB, MzCIPK, 
and MzGTL, the RT-PCR results of 14 tested genes (nearly 78%) had a good correlation with 
the microarray results. 
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We also selected ten genes not involved in transcriptional regulation for RT-PCR analysis. A gene of 
unknown function (Bankltl4981 14) and MzSDH (sorbitol dehydrogenase; Bankltl495577) were salt 
suppressed (Figure 2). The gene of unknown function went down sharply and reached the lowest 
expression at 1.5 h after treatment. In contrast to the previous reports, three genes (MzRCA; Ribulose-1, 
5-bisphosphate carboxylase; Bankltl495565), MzLhcb2 (Light-harvesting complex II protein; 
Bankltl495570), and MzRPE (Ribulose-phosphate 3-epimerase; Bankltl498105)) involved in 
photosynthesis increased their expression steadily after salt treatment. The expression of MzPIPl 
(Plasma Membrane Intrinsic Protein 1; Banklt 1494672), MzTIP2 (Tonoplast Intrinsic Protein 2; 
Bankltl494673), MzHSP (Heat Shock Protein 70; Bankltl498107), and MzRD22 (RD22; 
Bankltl495591) were also induced by salt stress. The expression peak of MzTIP2 appeared at eight 
hours, while the peak time of the other three genes appeared later than eight hours. A slight increase of 
expression was found for MzNTR under stress. 

2.2. Promoter Analysis 

The promoters of approximately 70 genes were unavailable (or having a promoter region <500 bp), as 
there are still gaps in apple genome as a result of many areas of repetitive DNA and high heterogeneity [11]. 
Of the 111 genes with promoters of over 2000 bp, salt-responsive elements were found that included 
ABRE (abscisic acid responsive element), ARE (anaerobic responsive element), TC-rich repeats 
(defense and stress responsive element), HSE (heat stress responsive element), DRE (drought 
responsive element), LTR (low-temperature responsive element), ERE (ethylene responsive element), 
and MBS (MYB bingding site). The average number of the cz's-elements per gene changed from 
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2.43 (ABRE) to 1 .0 (DRE). The top three cz's-elements with high frequency were ABRE, ARE, and MBS 
(Table 2). 



Table 2. The presentation of the cz's-elements in the salt-responsive gene in Malus zumi. 



Os-element 


ABRE 


ARE 


TC-rich 
repeats 


HSE 


DRE 


LTR 


ERE 


MBS 


Genes with the m-element 
among the Regulatory genes 


9/14 


13/14 


10/14 


9/14 


0/14 


8/14 


7/14 


10/14 


ROS elimination 


7/10 


5/10 


8/10 


8/10 


0/10 


4/10 


3/10 


7/10 
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Genes with 

Stress tolerance 

the 

Photosynthesis 

cM-element 

Transporter 

among the 

Metabolism 

functional 

Protein related 

genes 

Cell maintenance 
and development 


7/9 
9/14 


7/9 
10/14 


7/9 
10/14 


2/9 
6/14 


1/9 
2/14 


1/9 
2/14 


2/9 
2/14 


7/9 
12/14 


21/26 
7/9 
17/23 


20/26 
7/9 
19/23 


17/26 
6/9 
15/23 


13/26 
7/9 
16/23 


1/26 
0/9 
0/23 


12/26 
3/9 
10/23 


7/26 
1/9 
8/23 


18/26 
7/9 
18/23 


4/11 

2/9 


10/11 
7/9 


8/11 
4/9 


6/11 
7/9 


0/11 
0/9 


10/11 

3/9 


8/11 
1/9 


18/11 
2/9 


Average number of the 
cw-elements per gene 


2.43 


2.24 


1.75 


2.04 


1.00 


1.45 


1.19 


2.10 



2.3. Regulatory Network during Salt Stress 

In order to understand the molecular mechanism of apple salt response and explore the protein 
interaction network of Malus zumi under salt stress, an interaction network of salt responsive genes was 
constructed and the molecular regulatory pathways of apple were alluded to. A total of 155 unigenes, 
with annotated homologous genes in Arabidopsis, were included in an interaction map using 
Arabidopsis Interactions Viewer (http://bar.utoronto.ca (accessed on 22 June 201 1)) and visualized by 
Cytoscape 2.6.3 (http://www.cytoscape.org (accessed on 22 June 2011)). Some potential interactive 
proteins suggested by our salt stress responsive gene set were also included on the map. This suggested 
that the 58 genes presented in core set of interaction map played critical roles in salt stress response in 
Malus zumi, among which 28 genes were up-regulated and 13 genes were down-regulated from our data 
(Supplementary file (Table S2)). The interaction of the 58 proteins encoded by these genes 
(20 regulatory genes and 38 functional genes) will be affected by many aspects of plant activity that 
contribute together in salt resistance in Malus zumi (Figure 4). All of the three primary signaling 
pathways implicated in salt-stress response (CIPK/SOS, CDPK, and MAPK pathways) appear in apple, 
but it seems the CDPK pathway plays a central role in this response (Figure 3). 
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Figure 3. Interaction network of salt-responsive genes in Malus zumi. The interactions of 
salt responsive genes were predicted by their homologs in Arabidopsis. A total of 
22 up-regulated genes (red), 12 down-regulated genes (green) in microarray, and 15 
potentially interactive proteins (blue), were presented in the interaction network. CDPK 
appeared as quadrilateral, lines between CDPK and its interactive proteins were magenta; 
MAPK appeared in hexagon, lines between MAPK and its interactive proteins were khaki; 
CIPKs appeared in octagon, lines between CIPKs and its interactive proteins were red. The 
thick lines represented the direct interaction with kinases, the thin lines represented the 
interaction between proteins. Details of gene information were listed in Supplementary file 
(Table S2). 




3. Discussion 

3.1. Characteristics of Constructed cDNA Library 

A cDNA library was constructed from Malus zumi under salt stress by SMART ™ cDNA library 
construction kit. The titer of unamplified cDNA library was 2 x 10 6 pfu/mL. The recombination 
efficiency was approximately 90.5%. Colony PCR of 200 randomly selected clones showed that the 
insertions ranged from 250 bp to 2000 bp, and 115 insertions of different sizes were detected by 
polyacrylamide gel electrophoresis, which suggested the percentage of uniESTs of our cDNA library 
was around 58%. The percentage of the insertion sizes between 250 bp to 500 bp was approximately 3%; 
500 bp-750 bp, 28%; 750 bp-1000 bp, 36%, and 1.0-2.0 kb, 32%. To assess the percentage of the 
full-length cDNA insertions in our cDNA library, 30 recombinant clones were sequenced. A total of 
16 clones from the tested 30 clones (53%) had full-length cDNA sequences. The 5'UTR (untranslated 
regions) of full-length clones (0.7 kb-1 .6 kb) ranged from 70 bp to 300 bp and 3'UTR ranged from 1 50 bp to 
300 bp. Therefore, the summary statistics of the cDNA library suggested a high quality of cDNA. 
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3.2. Regulatory Network during Salt Stress 

An interaction network of salt responsive genes was constructed to describe the molecular 
mechanism of salt response in Malus zumi. Three primary signaling pathways in salt-stress response 
(CIPK/SOS, CDPK, and MAPK pathways) have been reported in plants [16]. In our network analysis, 
CDPK proteins appear to play central role in salt response in Malus zumi. They were tightly linked to 
proteins involved in osmolyte metabolism (sorbitol-dehydrogenase acting as key enzyme), and ROS 
scavenging [17], which are key strategies to defend against salt stress. CDPK was not found in our data, 
but HSP (Heat Shock Protein) was linked to CDPK and CIPK, central in the network, suggesting the 
possible role of CIPK in salt stress of Malus zumi. Similarly, MAPK were not identified in our study, but 
WRKY, Trihelix DNA-binding protein and AP2 appeared, which are potentially regulated by MAPKs. In 
summary, it seemed CDPK pathway functioned as the primary signaling pathway, CIPK and MAPK 
pathway might also be employed during salt stress response in Malus zumi. 

Although the sensor of salt stress in Malus zumi remains unknown, an osmosensor was deduced on 
the interaction map. This osmosensor was a histidine kinase and is highly homologous to a cytokinin 
receptor responsive to water deprivation in Arabidopsis [18]. As it is an aspect of salt stress, osmotic 
stress may be perceived by the osmosensor in Malus zumi to trigger the downstream molecular 
responses. Further efforts should be directed towards discovering the salt-stress sensors and identifying 
additional signaling components that mediate the salt-stress regulation of the expression and activities of 
ion transporters. 

3.3. Salt Stress Signaling 

Salt stress signaling played a central role in stress response. Kinases and phosphatases can mediate salt 
stress response through reversible phosphorylation of transcription factors and other functional genes [4]. A 
total of 16 genes involved in salt stress signaling were identified, seven genes were induced and nine genes 
were suppressed under salt stress (Table 3). There were 10 genes encoding transcription factors (zinc finger 
protein, IAA-LEUCINE RESISTANT3, WRKY, IAA26, Auxin response factor, GRAS family 
transcription factor, SALT TOLERANCE homolog protein, GT-like trihelix DNA-binding protein and two 
AP2 transcription factors) and 6 genes encoding kinases (two leucine-rich repeat protein kinase, tousled-like 
serine/threonine kinase, two CBL-interacting protein kinases and protein kinase family protein). 



Table 3. Signal elements involved in salt stress response. 

Number Expression a Putative annotation Genebank accession Identities p-value b 

Kinase 

1 I leucine-rich repeat transmembrane protein kinase NP_199948 64% 5.41 x 10 5 

1 S leucine-rich repeat family protein kinase NP_1 79336 42% 1.02 x 1(T 4 

1 I tousled-like serine/threonine kinase NP 568405 82% 3.20 x io~ 5 

1 I CIPK5 NP 568241 79% 7.74 x 10~ 5 

1 S CIPK6 NP 194825 80% 1.98 x 10~ 7 

1 S protein kinase family protein NP_194952 50% 2.25 x 1Q~ 5 
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Table 3. Cont. 

Number Expression a Putative annotation Genebank accession Identities p-value b 

Transcription factor 

1 I IAA-LEUCINE RESISTANT3 NP 200279 89% 4.97 x 10~ 4 

1 I IAA26 NPJ88271 80% 3.58 x 10~ 6 

1 S GT-like trihelix DNA-binding protein NP 177814 37% 3.04 x 10~ 6 

1 S zinc finger (CCCH-type) family protein NP 200670 59% 8.13 x 10~ 5 

1 I WRKY family transcription factor NP_001078015 50% 4.52 x 10~ 5 

1 S GRAS family transcription factor XP 002322514 52% 3.58 x 10~ 4 

2 S AP2 transcription factor NPJ73355 70% 4.26 x 10~ 6 
1 I SALT TOLERANCE homolog protein NP 849598 68% 5.37 x 10~ 7 
1 S Auxin response factor NP 182176 70% 2.78 x 1Q~ 3 



* S means suppression, I means induction; b p-value indicates probability of a gene showing significantly differentially expression between 
salt-treated samples and untreated samples at significance level of 0.05 using FDR correction. 

The WRKY, IAA-LEUCINE RESISTANT3, IAA26, and SALT TOLERANCE homolog protein 
genes were up-regulated by salt stress while five transcription factors were down-regulated under salt 
stress in the microarray data. This is the first report of WRKY responding to salt stress in apple. In 
contrast to a study in Arabidopsis [19], in which the non-induced AtSTO could enhance the salt 
tolerance of the transgenic plant by its over-expression, the expression of MzSTO increased under salt 
stress in our microarray result. This is also the first report of salt — inducible auxin responsive protein 
MzIAA and salt-repressed ARF (Auxin Response Factor) MzARF in apple. In Arabidipsis, IAA26 
belongs to Aux/IAA family, which binds ARF to inhibit its activation of downstream auxin response 
genes for plant growth [20,21]. Our results suggest that they may act collectively to retard growth under 
salt stress in Malus zumi. 

In the down-regulated transcription factors, the RT-PCR of MzSCL showed good correlation with the 
microarray result. The two AP2 transcription factors in our research belonged to DREB2 subfamily, 
which were expected to be salt stress inducible but was repressed in our microarray data. Our RT-PCR 
showed that they actually enhanced their expression under salt stress. 

3.4. Analysis of cis-Elements 

In the promoters of 1 1 1 salt responsive genes, ABRE motifs were present most frequently suggesting 
the major signaling pathways of salt response are probably ABA-dependent [22,23]. HSE were also 
enriched. The heat-shock protein appeared central of the interaction network, which interacted with HSF 
(binding to HSE) to release plants from abiotic stress [24,25]. HSP has an important role in salt response 
as well. MBS is a binding site for MYB transcription factors, which control many abiotic stress 
responses [26,27]. In addition, ARE motifs are involved in response to anaerobic conditions [28], LTR 
to low temperature and ERE to ethylene were present indicating the possible roles in low temperature, 
oxygen shortage, and ethylene response. 
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3.5. Cellular Response Strategies to Salt Stress in Malus zumi 

The transcriptional changes of salt-responsive genes allude to the signaling and cellular responsive 
strategies of salt stress in Malus zumi. A total of 235 salt-responsive unigenes were obtained whose 
expression changed more than four-fold after salt treatment. It has been reported in Arabidopsis, rice and 
other species that osmoprotection, ROS scavenging and ion homeostasis maintenance are the major 
cellular response strategies for plants under high-salinity environment [2,29,30]. Similar mechanisms 
were found in apple. Interestingly, possible new mechanism in Malus zumi was discovered, that involves 
sorbitiol. In the response of salt stress, osmoprotection plays a crucial role to maintain the cellular water 
potential and avoid water loss. A total of four genes with putative functions in osmoprotection were 
identified in Malus zumi. Sorbitol was reported to be an important osmolyte in Rosaceae plants [31]. 
Salt-inducible MzS6PDH (Sorbitol-6-Phosphate Dehydrogenase) and salt- repressed MzSDH were 
identified in our research. S6PDH is a well-known rate-limiting enzyme in sorbitol synthesis in mature 
leaves, while SDH converts sorbitol to fructose in sink organs [32]. It was conceivable that the 
down-regulation of MzSDH results in lower sorbitol consumption and up-regulation of MzS6PDH 
leading to more sorbitol synthesis and sorbitol accumulation under salt stress. Three sorbitol transporters 
of apple, which named MdSOT3, MdSOT4, and MdSOT5, were also induced by drought stress [33]. 
Mannitol was also a well-known osmolyte in many species [34,35]. Two genes encoding mannose 
6-phosphate reductase, the key enzyme in mannitol synthesis, were salt-inducible in our research. The 
accumulation of sorbitol and mannitol may be critical in osmoprotection under salt stress in Malus zumi. 

Under a high-salinity environment, water and iron transportation becomes important to controlling 
ion toxicity and maintaining water relations. An up-regulated MzNTR gene encoding ion channel 
involved in Na + /P symport, was identified. This transporter has been documented to enhance salinity 
tolerance in Arabidopsis [36]. A total of six aquaporins, transmembrane proteins facilitating water 
transportation across membranes, four plasma membrane intrinsic proteins and two tonoplast intrinsic 
proteins were obtained, among which five aquaporins were up-regulated by salt stress. MzPIPl;! was 
found to participate in salt stress response in Malus zumi [37]. 

Several ROS scavengers (glutathione transferase, peroxidase, catalase, and metallothionein) 
increased their expression. Interestingly an ABC transporter MzABC, potentially involved in phenolic 
transportation into the vacuole was salt responsive. Compared with model plant Arabidopsis, apple had 
large quantity of phenolic compounds. Some secondary metabolites including phenolic compounds 
were documented to be involved in ROS scavenging [38,39]. In addition, six proteins identified in our 
research were reported to be involved in ubiquitin-proteasome pathway in Arabidopsis [40-43], 
suggesting that ubiquitin-proteasome pathway for protein degradation is salt responsive. 

3.6. Stable Photosynthesis Suggests a Novel Salt Tolerance Mechanism in Malus zumi 

In other plant species high salinity inhibits photosynthesis [44-46]. However, in Malus zumi a total of 
35 genes from the identified 37 unigenes involved in photosynthesis increased their expression under 
salt stress (Table 4). All the genes encoding light harvesting complex proteins and carbonic anhydrase 
were induced in microarray result. Salt stress also enhanced the expression of MzRCA, which encoded a 
key enzyme in the Calvin cycle, suggesting this C0 2 -assimilating pathway was stimulated at high 



Int. J. Mol. Sci. 2013, 14 



21063 



salinity. In order to confirm the microarray results, genes functioning in light reaction {MzLhcbl), and 
dark reaction (MzRPE and MzRCA) were selected for semi-quantitative RT-PCR analysis, which 
showed consistent results with microarray data (Figure 3). Physiological research has shown stable 
photosynthesis in Malus zumi under salt stress [47]. Conceivably, the up-regulation of genes that play 
roles in light harvesting and CO2 uptake could result in stable photosynthesis under salt stress. Similar 
physiological responses of photosynthesis to high salinity have been documented in halotolerant plants, 
such as Styphnolobium japonicum and Atriplex centralasiatica [48,49]. The ability of Malus zumi to 
maintain normal photosynthetic activity suggests a possible novel salt tolerant mechanism. 



Table 4. Genes involved in photosynthesis in salt-stressed Malus zumi. 



Number 


Expression 


Putative annotation 


Genebank 
accession 


Identities 


p-value b 


1 


S 


1 -deoxy-D-xylulose-5-phosphate synthase 


NP_193291. 


41% 


5.20 x 


io- 4 


1 


I 


H13 protein 


NP 199522 


51% 


2.56 x 


10- 7 


2 


I 


photosystem II subunit R 


NPJ 78025 


75% 


3.43 x 


io- 6 


1 


I 


CHLOROPHYLL A/B BINDING PROTEIN 1 


NP 174286 


67% 


4.38 x 


IO" 3 


2 


I 


PSAF (photosystem I subunit F) 


NPJ 744 18 


72% 


3.75 x 


IO" 5 


1 


I 


photosystem II CP43 chlorophyll apoprotein 


YP 002149729 


90% 


6.54 x 


10~ 8 


1 


I 


photosystem I subunit D-2 


NP171812 


72% 


8.79 x 


io- 4 


1 


I 


Thylakoid membrane phosphoprotein of 14 kda 


NP_566086 


64% 


9.54 x 


IO" 7 


1 


I 


Photosystem I light harvesting complex gene 3 


NPJ 76347 


81% 


5.37 x 


IO" 6 


1 


I 


cytochrome b6 


NP_051088 


98% 


5.86 x 


10~ 7 


1 


I 


ATP synthase gamma chain 


NPJW265 


77% 


7.35 x 


IO" 4 


1 


I 


NADH dehydrogenase subunit 7 


NP_051115 


87% 


2.36 x 


10~ 7 


2 


I 


photosystem II 44 kDa protein 


NP_051055 


97% 


5.38 x 


IO" 5 


5 


I 


light-harvesting complex I protein Lhca3 


XP_002321218 


81% 


6.67 x 


io- 6 


1 


I 


Oxygen-evolving enhancer protein 


NP 201458 


83% 


4.31 x 


10~ 3 


1 


I 


PSBP-1 (PHOTOSYSTEM II SUBUNIT P-l) 


NPJ 72 153 


64% 


3.34 x 


IO" 4 


1 


I 


FERREDOXIN-NADP(+)-OXIDOREDUCTASE 1 


NP_201420 


76% 


4.53 x 


IO" 5 


1 




plastid-lipid-associated protein 


NPJ 923 11 


44% 


5.61 x 


IO" 4 


65 




ribulose- 1 ,5-bisphosphate carboxylase 


CAA79857 


91% 


5.64 x 


io- 6 


1 




photosystem I P700 chlorophyll a apoprotein A2 


NP 051058 


98% 


6.35 x 


io- 3 


7 




light-harvesting complex II protein Lhcb2 


XP_002321186 


92% 


6.89 x 


IO" 7 


1 




light-harvesting complex I protein Lhca4 


XP_002330127 


75% 


8.32 x 


io- 6 


1 




light-harvesting complex I protein Lhca2 


XP 002299309 


94% 


3.57 x 


IO" 4 


1 




PHOTOSYSTEM II SUBUNIT 0-2 


NPJ 90651 


75% 


4.25 x 


io- 5 


1 




PsbP domain-containing protein 


NP 565131 


65% 


6.31 x 


IO" 7 


5 




light-harvesting complex II protein Lhcb 1 


XP 0023 16737 


95% 


2.58 x 


io- 2 


6 




light-harvesting complex II protein Lhcb6 


XPJ02303160 


82% 


5.34 x 


io- 5 


1 




light-harvesting complex II protein Lhcb5 


XP 002329192 


81% 


6.41 x 


IO" 7 


2 




photosystem II 22 kDa protein 


NP 00 1150026 


91% 


3.59 x 


io- 3 


2 




protochlorophyllide reductase 


NP 200230 


65% 


5.34 x 


io- 7 


1 




Clp protease proteolytic subunit 6 


NP_563893 


45% 


4.37 x 


io- 5 


2 




magnesium chelatase H subunit 


AC057443 


100% 


1.82 x 


io- 8 


1 




Photosystem I P700 chlorophyll A apoprotein 


NP 00 1044491 


60% 


2.59 x 


io- 4 


1 




carbonic anhydrase 


NPJ 86799 


61% 


1.56 x 


io- 4 


1 




starch synthase 


NP 174566 


36% 


8.24 x 


IO" 5 


1 




glucose- 1-phosphate adenylyltransferase 


NPJ 97423 


81% 


2.34 x 


io- 7 


1 




ribulose-phosphate 3-epimerase,RPE 


NP 200949. 


83% 


5.20 x 


IO" 4 



a S means suppression, I means induction; /?-value indicates probability of a gene showing significantly 
differentially expression between salt-treated samples and untreated samples at significance level of 0.05 using 
FDR correction. 
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3. 7. Epigenetic and Novel Regulation in Salt Stress Response 

Histone modifications play a significant role in stress response [50]. The genes encoding histone 
deacetylase and histone-lysine iV-methyltransferase were strongly repressed in our microarray results, 
suggesting involvement in salt response. The putative interaction between the two proteins and other 
salt-responsive proteins were presented in the protein interaction network. 

About 21% (50 unigenes) of salt responsive genes identified in our study were functionally unknown. 
A total of 24 novel sequences, showing no similarity to any sequence in public databases, were 
identified. The expression patterns of one functionally unknown gene (MzUKl) were confirmed with 
RT-PCR. These genes, with no known matches in other species or homologous to other genes with 
known function, may be useful for gene mining to detect unique salt responsive pathways of other 
woody species. 

A total of 16 unigenes involved in other biotic and abiotic stresses were presented in our results, 
among which five genes were homologous to early drought induced protein, disease resistance protein, 
pathogenesis-related thaumatin family protein, heat shock protein 70 and UV-induced protein [51-53] 
(Supplementary file (Table SI)). This suggests possible cross-talk between salt stress and other stresses. 

4. Experimental Section 

4.1. Plant Materials 

All the seedlings of Malus zumi were micropropagated in vitro by culturing in MS medium 
supplemented with 0.5 mg/L BA. Seedlings were transferred to rooting medium (1/2 MS + 1.0 mg/L 
IAA) when they were at least 3 cm tall. Rooted plants were transplanted into in 1/2 nutrient solution. 
Nearly one month later, new roots grew out and the rooted plantlets were transferred to 
complete nutrient solution in a growth chamber (800 jamol ■ m • s light intensity with 1 2 h photoperiod) 
as described by Han et al. (1994). The day/night temperature and humidity were maintained at 22-25 °C/ 
15-18 °C and 45~50%/60~70%, respectively. Seedlings higher than 10 cm were transferred to nutrient 
solution supplemented with 150 mM NaCl for salt treatment, nutrient solution without NaCl was used 
for control. Two biological replicate samples of the leaves and roots were collected at 0 h, 
0.5 h, 1 h, 1.5 h, 3 h, 6 h, 12 h, 1 day, 3 day, 5 day, and 7 day after salt treatment, and each sample 
contained five plants. To prepare RT-PCR materials, Seedlings higher than 10 cm were transferred to 
nutrient solution supplemented with 150 mM NaCl for salt treatment. Leaves and roots were collected at 
0 h, 2 h, 4 h, and 8 h. These materials were stored at -80 °C until use. 

4.2. RNA Extraction and Construction of the Normalized Full-Length cDNA Library o/Malus zumi 

Total RNA was extracted separately from stored samples, using the CTAB method [54]. Equal 
amount of total RNAs from salt- treated leaves and roots (0 h, 0.5 h, 1 h, 1 .5 h, 3 h, 6 h, 12 h, 1 day, 3 day, 
5 day, and 7 day) were mixed to isolate mRNA by Poly Attract mRNA isolation system III (Promega) for 
cDNA library construction. The leaves and roots from control samples were collected and mixed at the 
same time point as treated samples. cDNA was purified and normalized to reduce the redundancy. 
The cDNA library was constructed by using CLONTECH SMART ™ cDNA Library Construction Kit. 
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The primary cDNA library was titered and amplified following the protocols. To check the 
recombination efficiency and insertion size, 200 colony PCRs were conducted to amplify the inserts 
using the commercial primers (Ml 3-48). 

4.3. Preparation of Slides for cDNA Microarrays 

A total of 15,000 cDNA clones were randomly selected in the cDNA microarray analysis. Insertions 
of the cDNA clones were amplified by PCR using commercial primers (M13-48). The yield and 
amplification quality of PCR products were checked on 1% (w/v) agarose gel. Purified PCR fragments 
were arrayed from 3 84- well micro titer plates onto a poly-l-Lys-coated micro slide glass by 
Cartesian7500 spotting robotics (Cartesian Inc., Newton, MA, USA). Each fragment was dotted twice 
on the same filter. After spotting, the slides were hydrated for two hours and dried for 0.5 h, which was 
followed by UV cross-linking and treatment with 0.2% SDS, H 2 0 and 0.2% NaBFL. The slides were 
dried again before use. 

Salt-responsive cDNA candidates were screened based on these criteria: (1) signal intensity changed 
more than two-fold on average in both dye-swap experiments; (2) /»-values of cDNA candidates were 
less that 0.05; and (3) the spot did not give irregular signals such as due to deformation or dust. 

4.4. Probe Labeling, Slide Hybridization and Scanning 

Two biological replicate samples of the leaves and roots were collected at Oh, 0.5 h, 1 h, 1.5 h, 3 h, 
6 h, 12 h, 1 day, 3 day, 5 day, and 7 day after salt treatment. Each sample represented five plants for RNA 
extraction. The control samples were always same total RNA mixture from five untreated plants. Equal 
amount of Total RNA extracted from salt-stressed samples were mixed together and labeled with Cy5 by 
reverse transcription. The cDNA derived from untreated control sample, which was labeled with Cy3 by 
reverse transcription, was used as the expression reference. The reverse transcription cocktail included 
200 U Superscript RNase H" reverse transcriptase (Gibco BRL, Rockville, MD, USA). The Cy3- or 
Cy5-labeled cDNAs were dissolved in 20 jxL solution, containing 5x SSx and 0.2% SDS. After 
denaturation at 95 °C, chips were placed in ethanol for 30 s. The probes were placed on ice respectively. 
Hybridizations were performed for 16-18 h at 42 °C in Hybricasette (Shellab Ltd., Cornelius, OR, 
USA). The slides were then transferred into 0.1% SSC and shaken gently for 20 min. After the slides 
were rinsed twice, they were spun in a centrifuge and dried by blowing. Dye swap were used in all 
experiments to compare gene expression changes. The fluorescent signatures were captured using a 
ScanArray 3000 (GSI; Luminomics, Billercia, MA, USA). 

4.5. Data Analysis of Microarray Results 

Image analysis and signal quantification were performed using ImaGene III software (BioDiscovery, 
Los Angeles, CA, USA). Background subtraction was carried out using the average of the lowest 5% of 
spot fluorescence intensities. The log2 ratios (signal intensity of Cy5/Cy3) were normalized using the 
intensity-based Loess method with R language. Wilcoxon Signed Rank Test was applied to compute p 
value using multtest packages in bioconductor (http://www.bioconductor.Org/packages/2. 12/bioc/html/ 
multtest.html) [55-57]. The False Discovery Rate (FDR) was used for multiple comparison corrections 
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according to Benjamini-Hochberg method [58]. The differentially expressed genes were defined by a 
log-scale ratio with a corrected p value < 0.05. The ratio of intensities between treatment and control 
indicated gene transcription. When the change of signal intensity between treatment and control 
exceeded four-fold, the corresponding genes were defined as salt-responsive genes. Then the 
salt-responsive genes were sequenced (BGI inc. Beishan Street Beishan industrial estate, Yan Tian 
district, Shenzhen, China) and their homologies were screened with in the GenBank database using 
BLASTx search program (http://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 5 April 2010)). The 
unigenes were grouped into functional categories according to Gene Ontology (GO) prediction. 

4.6. RT-PCR Analysis 

Reverse transcription was carried out with total RNA (1 (j,g) isolated from stressed sample/control 
using SMART™ PCR cDNA Synthesis kit following manufacturer's instructions. The primers for each 
selected gene were shown in Supplementary file (Table S3). PCRs were performed for one cycle at 
94 °C for 3 min, followed by 28 to 38 cycles of 94 °C for 30 s, 52 °C to 53 °C for 30 s, and 72 °C for 40 
s, with a final extension of 10 min at 72 °C. Mz^C77jV(AB638619.1) was used as the internal control. 
The relative amounts of PCR products were quantified by direct scanning of ethidium bromide-stained 
1 % TAE-agarose gels with Alpha Imager imaging system equipped with the AlphaEaseFC Windows 
Software. The expression levels of tested genes were normalized according to the corresponding 
MzACTIN amplifications. 

4. 7. Promoter Analysis 

A total of 1 85 arabidopsis genes homologous to annotated apple genes were examined using the GDR 
database (http://www.rosaceae.org/gb/gbrowse/malus_x_domestica/ (accessed on 10 April 2010)), 
allowing the promoter sequence of 1 1 1 genes to be obtained. The 2000bp fragment upstream of the start 
codon was assumed to be promoter region. The czs-elements of promoter sequences were analyzed by 
PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ (accessed on 12 April 2010)). 

4.8. Exploring the Interaction Network 

Unigenes identified as responding to salt stress in the microarray experiment were searched against 
the TAIR database for homologous genes in Arabidopsis. Interactions of these genes were predicted by 
Arabidopsis Interactions Viewer (http://bar.utoronto.ca) and visualized by Cytoscape 2.6.3 
(http://www.cytoscape.org (accessed on 22 June 2011)). The predicted interacting proteins that might 
play roles in plant salt stress response were then presented in the network. 

5. Conclusions 

We report a global expression profile of salt-responsive genes by microarray analysis in apple 
rootstock Malus zumi. In contrast to previous reports in other species, we found increased expression of 
genes involved in photosynthesis under salt stress and new mechanisms for ROS scavenging and 
osmoprotection. An interaction framework of salt stress responsive genes was generated to summarize 
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the salt stress response of Malus. Our work contributes to furthering the understanding of salt response 
mechanisms in apple trees and engineering apple plants with enhanced salt tolerance. 

Acknowledgments 

Our research was supported by the 973 Program (Project No.: 2012CB1 14200), Beijing Natural 
Science Foundation (Project No. :61 12012). 

Conflicts of Interest 

The authors declare no conflict of interest. 
References 

1. Munns, R.; Tester, M. Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 2008, 59, 651-681. 

2. Zhu, J.K. Regulation of ion homeostasis under salt stress. Curr. Opin. Plant Biol. 2003, 6, 441^145. 

3. Hasegawa, P.M.; Bressan, R.A.; Zhu, J.K.; Bohnert, H.J. Plant Cellular and Molecular Responses 
to High Salinity. Rev. Plant Physiol. Plant Mol. Biol. 2000, 51, 463-499. 

4. Zhu, J.K. Salt and drought stress signal transduction in plants. Annu. Rev. Plant Biol. 2002, 53, 
247-273. 

5. Hu, H.; Dai, M.; Yao, J.; Xiao, B.; Li, X.; Zhang, Q.; Xiong, L. Overexpressing aNAM, ATAF, and 
CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc. Natl. 
Acad. Sci. USA 2006, 103, 12987-12992. 

6. Shi, H.; Lee, B.H.; Wu, S.J.; Zhu, J.K. Overexpression of a plasma membrane Na+/H+ antiporter 
gene improves salt tolerance in Arabidopsis thaliana. Nat. Biotechnol. 2003, 21, 81-85. 

7. Wen, X.P.; Pang, X.M.; Matsuda, N.; Kita, M.; Inoue, H.; Hao, Y.J.; Honda, C; Moriguchi, T. 
Over-expression of the apple spermidine synthase gene in pear confers multiple abiotic stress 
tolerance by altering polyamine titers. Transgenic Res. 2008, 17, 251-263. 

8. Yin, R.; Bai, T.H.; Ma, F.W.; Wang, X.J.; Li, Y.H.; Yue, Z.Y. Physiological responses and relative 
tolerance by Chinese apple rootstocks to NaCl stress. Sci. Hortic. 2010, 126, 247-252. 

9. Flowers, T.J. Salinisation and horticultural production. Sci. Hortic. 1999, 78, 1-4. 

10. Xia, Y.; Huimin, H.; Shu, H.R.; Wang, T.M.; Liu, D.X.; Fang, Y.F.; Chai, C.H. Changes of leaf 
membrane penetration, proline and mineral nutrient contents of young apple tree under NaCl stress. 
J. Fruit Sci. 2005, 22, 1-5. 

11. Velasco, R.; Zharkikh, A.; Affourtit, J.; Dhingra, A.; Cestaro, A.; Kalyanaraman, A.; Fontana, P.; 
Bhatnagar, S.K.; Troggio, M.; Pruss, D.; et al. The genome of the domesticated apple 
(Malus x domestica Borkh.). Nat. Genet. 2010, 42, 833-839. 

12. Janssen, B.J.; Thodey, K; Schaffer, R.J.; Alba, R.; Balakrishnan, L.; Bishop, R.; Bowen, J.H.; 
Crowhurst, R.N.; Gleave, A.P.; Ledger, S.; et al. Global gene expression analysis of apple fruit 
development from the floral bud to ripe fruit. BMC Plant Biol. 2008, 8, 16. 

13. Kawaura, K; Mochida, K; Ogihara, Y. Genome-wide analysis for identification of salt-responsive 
genes in common wheat. Funct. Integr. Genomics 2008, 8, 277-286. 



Int. J. Mol. Sci. 2013, 14 



21068 



14. Zhao, J.; Wang, J.; An, L.; Doerge, R.W.; Chen, Z.J.; Grau, C.R.; Meng, J.; Osborn, T.C. Analysis 
of gene expression profiles in response to Sclerotinia sclerotiorum in Brassica napus. Planta 2007, 
227, 13-24. 

15. Wei, G.; Tao, Y.; Liu, G.; Chen, C; Luo, R.; Xia, H.; Gan, Q.; Zeng, H.; Lu, Z.; Han, Y.; et al. 
A transcriptomic analysis of superhybrid rice LYP9 and its parents. Proc. Natl. Acad. Sci. USA 
2009,106, 7695-7701. 

16. Xiong, L.; Schumaker, K.S.; Zhu, J.K. Cell signaling during cold, drought, and salt stress. 
Plant Cell 2002, 14, S165-S183. 

17. Battaglia, C; Hsu, CM.; Soderstrom, K.; Escarre, J.; Haug, F.J.; Charriere, M.; Boccard, M.; 
Despeisse, M.; Alexander, D.T.; Cantoni, M.; et al. Light trapping in solar cells: Can periodic beat 
random? ACS Nano 2012, 6, 2790-2797. 

18. Tran, L.S.; Urao, T.; Qin, F.; Maruyama, K.; Kakimoto, T.; Shinozaki, K.; Yamaguchi-Shinozaki, K. 
Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to 
abscisic acid, drought, and salt stress in Arabidopsis. Proc. Natl. Acad. Sci. USA 2007, 104, 
20623-20628. 

19. Nagaoka, S.; Takano, T. Salt tolerance-related protein STO binds to a Myb transcription factor 
homologue and confers salt tolerance in Arabidopsis. J. Exp. Bot. 2003, 54, 2231-2237. 

20. Ding, D.; Zhang, L.F.; Wang, H.; Liu, Z.J.; Zhang, Z.X.; Zheng, Y.L. Differential expression of 
miRNAs in response to salt stress in maize roots. Ann. Bot. 2009, 103, 29-38. 

21. Teale, W.D.; Paponov, I. A.; Palme, K. Auxin in action: Signalling, transport and the control of 
plant growth and development. Nat. Rev. Mol. Cell. Biol. 2006, 7, 847-859. 

22. Mundy, J.; Yamaguchi-Shinozaki, K.; Chua, N.H. Nuclear proteins bind conserved elements in the 
abscisic acid-responsive promoter of a rice rab gene. Proc. Natl. Acad. Sci. USA 1990, 87, 
1406-1410. 

23. Ross, C; Shen, Q.J. Computational prediction and experimental verification of HVAl-like abscisic 
acid responsive promoters in rice (Oryza sativa). Plant Mol. Biol. 2006, 62, 233-246. 

24. von Koskull-Doring, P.; Scharf, K.D.; Nover, L. The diversity of plant heat stress transcription 
factors. Trends Plant Sci. 2007, 12, 452-457. 

25. Wang, W.; Vinocur, B.; Shoseyov, O.; Altman, A. Role of plant heat-shock proteins and molecular 
chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9, 244-252. 

26. Dai, X.; Xu, Y.; Ma, Q.; Xu, W.; Wang, T.; Xue, Y.; Chong, K. Overexpression of an R1R2R3 
MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic 
Arabidopsis. Plant Physiol. 2007, 143, 1739-1751. 

27. Martin, C; Paz-Ares, J. MYB transcription factors in plants. Trends Genet. 1997, 13, 67-73. 

28. Dolferus, R.; Jacobs, M.; Peacock, W.J.; Dennis, E.S. Differential interactions of promoter 
elements in stress responses of the Arabidopsis Adh gene. Plant Physiol. 1994, 105, 1075-1087. 

29. Ashraf M, F.M. Roles of glycine betaine and proline in improving plant abiotic stress resistance. 
Environ. Exp. Bot. 2007, 59, 206-216. 

30. Tsugane, K.; Kobayashi, K.; Niwa, Y.; Ohba, Y.; Wada, K.; Kobayashi, H. A recessive arabidopsis 
mutant that grows photoautotrophically under salt stress shows enhanced active oxygen 
detoxification. Plant Cell 1999, 11, 1195-1206. 



Int. J. Mol. Sci. 2013, 14 



21069 



31. Kanayama, Y.; Moriguchi, R.; Yamaki, S.; Deguchi, M.; Kanahama, K. Effects of environmental 
stresses and abscisic acid on Sorbitol-6-Phosphate dehydrogenase expression in Rosaceae fruit 
trees. Acta Hortic 2007, 375-381. 

32. Yamaki, S.I.K. Role of four sorbitol related enzymes and invertase in the seasonal alteration of 
sugar metabolism in apple tissue. J. Am. Soc. Hort. Sci. 1986, 111, 134-137. 

33. Li, F.; Lei, H.J.; Zhao, X.J.; Tian, R.R.; Li, T.H. Characterization of Three Sorbitol Transporter 
Genes in Micropropagated Apple Plants Grown under Drought Stress. Plant Mol. Biol. Rep. 2012, 
30, 123-130. 

34. Pharr, D.M.; Stoop, J.M.H.; Williamson, J.D.; Feusi, M.E.S.; Massel, M.O.; Conkling, M.A. The 
Dual Role of Mannitol as Osmoprotectant and Photoassimilate in Celery. Hortscience 1995, 30, 
1182-1188. 

35. Tarczynski, M.C.; Jensen, R.G.; Bohnert, H.J. Stress protection of transgenic tobacco by 
production of the osmolyte mannitol. Science 1993, 259, 508-510. 

36. Cubero, B.; Nakagawa, Y.; Jiang, X.Y.; Miura, K.J.; Li, F.; Raghothama, K.G.; Bressan, R.A.; 
Hasegawa, P.M.; Pardo, J.M. The phosphate transporter PHT4;6 is a determinant of salt tolerance 
that is localized to the Golgi apparatus of Arabidopsis. Mol. Plant 2009, 2, 535-552. 

37. Li, Q.T.; Zhao, Y.; Zhou, Q.Y.; Huang, L.F.; Jiang, Y.Z.; Chen, H.; Kong, J. Cloning and 
expression analysis of MzPIPl;l gene from Malus zumi Mats. J. China Agric. Univ 2012, 17, 6. 

38. Bartholomew, D.M.; Van Dyk, D.E.; Lau, S.M.; O'Keefe, D.P.; Rea, PA.; Viitanen, P.V. Alternate 
energy-dependent pathways for the vacuolar uptake of glucose and glutathione conjugates. Plant. 
Physiol. 2002, 130, 1562-1572. 

39. Marrs, K.A.; Alfenito, M.R.; Lloyd, A.M.; Walbot, V. A glutathione S-transferase involved in 
vacuolar transfer encoded by the maize gene Bronze-2. Nature 1995, 375, 397-400. 

40. Balakirev, MY.; Tcherniuk, S.O.; Jaquinod, M.; Chroboczek, J. Otubains: A new family of 
cysteine proteases in the ubiquitin pathway. EMBO Rep. 2003, 4, 517-522. 

41. Hanna, J.; Meides, A.; Zhang, D.P.; Finley, D. A ubiquitin stress response induces altered 
proteasome composition. Cell 2007, 129, 747-759. 

42. Jentsch, S. The ubiquitin-conjugation system. Annu. Rev. Genet. 1992, 26, 179-207. 

43. Singer, T.; Haefner, S.; Hoffmann, M.; Fischer, M.; Ilyina, J.; Hilt, W. Sit4 phosphatase is 
functionally linked to the ubiquitin-proteasome system. Genetics 2003, 164, 1305-1321. 

44. Chaves, M.M.; Flexas, J.; Pinheiro, C. Photosynthesis under drought and salt stress: Regulation 
mechanisms from whole plant to cell. Ann. Bot. 2009, 103, 551-560. 

45. Dionisio-Sese, M.L.; Tobita, S. Effects of salinity on sodium content and photosynthetic responses 
of rice seedlings differing in salt tolerance. J. Plant Physiol. 2000, 157, 54-58. 

46. Downton, W.J.S. Photosynthesis in Salt-Stressed Grapevines. Aust. J. Plant. Physiol. 1977, 4, 
183-192. 

47. Ma, L.Q.; Han, Z.H.; Zhou, E.F.; Xu, X.F. Effect of salt stress on protective enzyme system of 
membrane in Malus zumi and M.baccata. J. Fruit Sci. 2006, 23, 495-499. 

48. Qiu, N.W.; Lu, Q.T.; Lu, CM. Photosynthesis, photosystem II efficiency and the xanthophyll cycle 
in the salt-adapted halophyte Atriplex centralasiatica. New Phytol. 2003, 159, 479-486. 

49. Zhang, C.H.; Shen,Y.B.; Yin, W.L.; Pan, Q.H.; Zhao, Y.G. Effect of salt stress on photosynthesis 
and growth of four tree species seedings. Sci. Silvae Sin. 2002, 38, 27-31. 



Int. J. Mol. Sci. 2013, 14 



21070 



50. Chinnusamy, V.; Zhu, J.K. Epigenetic regulation of stress responses in plants. Curr. Opin. Plant 
Biol. 2009, 12, 133-139. 

51. Davey, S.; Nass, M.L.; Ferrer, J.V.; Sidik, K.; Eisenberger, A.; Mitchell, D.L.; Freyer, G.A. The 
fission yeast UVDR DNA repair pathway is inducible. Nucleic Acids Res. 1997, 25, 1002-1008. 

52. Lee, G.J.; Vierling, E. A small heat shock protein cooperates with heat shock protein 70 systems to 
reactivate a heat-denatured protein. Plant Physiol. 2000, 122, 189-198. 

53. Stintzi, A.; Heitz, T.; Prasad, V.; Wiedemannmerdinoglu, S.; Kauffmann, S.; Geoffroy, P.; Legrand, 
M.; Fritig, B. Plant Pathogenesis-Related Proteins and Their Role in Defense against Pathogens. 
Biochimie 1993, 75, 687-706. 

54. Zhang, Y.G.; Kong, J.; Wang, Y.; Xu, X.F.; Liu, L.L.; Li, T.Z.; Zhu, Y.J.; Han, Z.H. Isolation and 
characterisation of a nicotianamine synthase gene MxNasl in Malus xiaojinensis. J. Hortic. Sci. 
Biotech. 2009, 84, 47-52. 

55. Wilson, I.W.; Kennedy, G.C.; Peacock, J.W.; Dennis, E.S. Microarray analysis reveals vegetative 
molecular phenotypes of Arabidopsis flowering-time mutants. Plant Cell Physiol. 2005, 46, 
1190-1201. 

56. Mohr, P.G.; Cahill, D.M. Suppression by ABA of salicylic acid and lignin accumulation and the 
expression of multiple genes, in Arabidopsis infected with Pseudomonas syringae pv. tomato. 
Funct. Integr. Genomics 2007, 7, 181—191. 

57. Pollard, K.S.; Dudoit, S.; van der Laan, M.J. Multiple Testing Procedures: R Multtest Package and 
Applications to Genomics. Bioinformatics and Computational Biology Solutions Using R and 
Bioconductor; Springer: Berlin, Germany, 2005; Chapter 15. 

58. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate — A Practical and Powerful 
Approach to Multiple Testing. J. Roy. Stat. Soc. B 1995, 57, 289-300. 



© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



